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The influence of surface anisotropy on the magnetization processes of maghemite nanoparticles with ellip- 
soidal shape is studied by means of Monte Carlo simulations. Radial surface anisotropy is found to favor 
the formation of hedgehog-like spin structures that become more stable as the surface anisotropy constant 
at the surface ks is increased form the value at the core. We have studied the change in the low temperature 
hysteresis loops with the particle aspect ratio and with ks, finding a change in the magnetization reversal 
mode as ks or the particle elongation is increased. 
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1 Introduction 

Most spins in a magnetic nanoparticle are at the particle surface where the reduced coordination and the 
broken translational symmetry of the lattice are responsible for changes in some local physical properties 
of the material. Among them, there is the reduction in exchange interaction constants and the increased 
local anisotropy at the particle surface with respect to bulk values 1 1 1 phenomenological model that de- 
scribed this last phenomenon was first introduced by Neel in \ 2\, but experimental determination of the 
microscopic structure and properties of the surface of nanoparticles is still a challenge to be solved. For 
example, an estimate of the surface anisotropy constant can only extracted indirectly from magnetic mea- 
surements by fitting to an effective anisotropy constant |3|. However, there are increasing evidences from 
experiments that surface effects play a dominant role in the anomalous low temperature magnetic prop- 
erties of nanoparticle systems |j4i. Therefore, simulations of magnetic properties in which microscopic 
parameters extracted form experiments are used |3[6|[7] mav serve as a benchmark to test the phenomeno- 
logical models at hand. In this article, we use the Monte Carlo method to simulate the magnetization 
processes of individual nanoparticles modeled by classical Heisenberg spins on the maghemite lattice with 
radial surface anisotropy kg different from that on the core. 

2 Model and simulation method 

The model considered is an extension to Heisenberg spins of our previous MC simulation for an Ising spin 
model of a 7-Fe2C<3 (maghemite) particle 1 8 1. The particles considered have spherical or ellipsoidal shape 
with diameter D = 3 and elongations L = 3 — 8 (in units of the cell size a), containing from TV = 347 
to 883 spins. The surface spins of the particle are considered to be those on the outer unit cell and have 
anisotropy constant ks and radial local anisotropy axes along the directions rii, whereas core spins have 
uniaxial anisotropy kc pointing along the z-axis. Indirect estimations of the anisotropy constants yield 
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values in the range kc = 0.01 — 1 K, ks = 1 — 100 K, so we will fix kc — 1 K and vary only ks- The 
interaction Hamiltonian in temperature units can be written as 

H/k B = -J2 J ijS i -S j -J2h-S i - £ [kc{Sl) 2 ]+ J2 [ fc s(^-n,) 2 ] ■ (1) 

(hi) i=1 i=1 



In the first term, Jy are the nearest-neighbours exchange constants among the spins in the two sublattices 
with different coordination of maghemite, all of them are antiferromagnetic, being Jto the strongest in 
magnitude and making maghemite a ferrimagnetic compound |9|. We have computed equilibrium config- 
urations by starting from a disordered state at a high temperature and decreasing the temperature in steps 
dT — 1. Hysteresis loops are started from a saturating field along the core easy-axis that is decreased 
in steps dh = —1, during which the magnetization is averaged over 1000 MC steps. In both cases, the 
standard MC method for continuous spins has been used with a trial jump consisting in a combination of 
uniform and small angle around the current spin direction that allows to obtain reliable simulations for treat 
the weak and strong anisotropy cases with the same algorithm, see [ 1 1 for the details. 

Fig. 1 Equilibrium spin configurations of an ellipsoidal maghemite particle with radial surface anisotropy, diameter 
D — 3a and elongation L = 8a, achieved after cooling form high temperature to T = 0. Upper (lower) panels are 
cross-sections along the equatorial plane (long-axis) showing only the spins on the central unit cell. 



3 Results 

In Fig. n we show snapshots of the spin configurations attained after the cooling process previously 
described, for several values of the surface anisotropy constant ks. The pictures present a cut of one unit 
cell width along the equatorial plane and along the major axis of the particle. For values of ks smaller 
than a certain critical value k s (which depends on the particle size and on the shape), the radial anisotropy 
at the surface is not strong enough to cant the core spins towards the radial direction. Therefore, the core 
spins are antiferromagnetically aligned with a small net magnetization component along the z axis due a 
small number of non-compensated spins (see Fig. \J]for ks = 10). However, at the particle surface, the 
spins present already a high degree of disorder due to the competition between the reduced local exchange 
field (locally reduced with respect to that on the core spins) and the anisotropy, inducing some degree of 
frustration at the surface. As ks is further increased (see Fig. [^for ks = 20), the disorder of the surface 
spins is transmitted to the particle core and results in the departure of core spins from the local anisotropy 
axis towards the radial direction. Finally, when anisotropy dominates over exchange energy (ks > kg, 
see Fig. [^for ks — 100), the hedgehog-like structures are formed, with the spins pointing along the local 
radial direction and having antiparallel orientations with the nearest neighbors on the other sublattice. 

We have also studied the influence of surface anisotropy and particle shape on the magnetization reversal 
processes. For this purpose, zero temperature hysteresis loops have been simulated using the procedure 
described previously for ellipsoidal particles with elongations L = 3 — 8 and diameter D = 3. Some 
selected results are presented in Fig. [2] where the total magnetization component along the field direction 
(left panels) is shown together with the contribution of the surface and core spins to M n = J2iLi &i ' ^» 
, the sum of the projections of spins onto the local anisotropy axis (central and right panels). By comparing 
the loops for ks= 10 to those for ks= 100, we see that the coercive field increases and the remanence 
decreases with increasing surface anisotropy, independently of the elongation L of the particle. Moreover, 
the presence of disordered groups of spins at the surface induced by surface anisotropy makes the loops 
more elongated and increases the closure fields of the loops as found also in experiments on ferrimagnetic 
oxides II II . The rounding of the loops near the coercive field clearly indicates a progressive departure from 
a uniform reversal mechanism with increasing ks. When looking only at the M z component, not much 
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Fig. 2 Hysteresis loops for ellipsoidal particles with two values of the radial surface anisotropy constant ks = 10, 100. 
Dashed lines are for a spherical particle with diameter D = 3, red circles and blue squares for an ellipsoidal particle 
with D = 3 and L — 6, 8 respectively. Left panels display the total magnetization while central and right panels 
show the surface and core contributions to M n , the sum of projections of the magnetization in to the local easy-axis 
directions. 



difference is appreciated between the loops for particles with different L because of the compensation of 
the spin components transverse to the field direction due to the symmetry of revolution of the particles 
around the z axis. However, upon further inspection of the M^ ur ^ and M^ ore contributions and animated 
snapshots taken along the loops 1 12 1, the details of the reversal process can be better under-stood. Let 
us notice first that, when ks is increased from k s , the reversal mechanism changes from quasi-uniform 
(induced by the core) rotation at low ks values, to a process in which the formation of surface hedgehog- 
like structures induce the non-uniform switching of the whole particle. In the first regime (ks = 10 case 
in Fig. |3, the core and surface spins point mostly along the z axis (M^ ore ~ 1, M^ ur f <^ 1) except 
near the coercive field, where they make short excursions towards the radial direction driven by the surface 
anisotropy (see the dips in M^ ore and the cusps in M^ ur ^). However, for ksi k s (ks = 100 case in Fig. 
[3, surface spins remain near the local radial easy-directions (M^ ur ^ ~ 1) during the reversal, while the 
core spins are dragged out the z local easy-axis by the surface spins during the reversal, except for values 
of h near the closure field. This is indicated by the widening of the dips in M„ re and the global decrease 
of M^ ore values as ks increases. Finally, let us remark also that, for all the ks considered, the M^ ur ^ 
values along the whole hysteresis loops increase with increasing L, which indicates that the surface spins 
remain closer to the local radial direction during the reversal as particles become more elongated. When 
increasing L, the dips in M^ ore become less profound for ks > k s , and indication that reversal of core 
spins along the radial direction is suppressed by the elongation. However, for weak anisotropy (ks < kg), 
the more elongated the particles are, the greater the deviation of surface spins towards the radial direction 
during the reversal. 
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